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We consider the collective mode spectrum of a normal Fermi gas in a spherical harmonic trap. 
Using a self-consistent random-phase-approximation, we systematically examine the effects of the 
two-body interactions on the modes of various symmetries. For weak coupling where the spectrum 
is shifted very little away from the ideal gas case, a sum-rule approach is shown to work well. For 
stronger coupling, the interplay between the single particle and the collective excitations causes 
effects such as mode splitting and Landau damping. A finite low frequency response present at 
stronger coupling is predicted for the quadrupole mode. Finally, we briefly discuss the effect of a 
finite temperature on the spectrum and the excitation of higher collective modes. 
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^""^ , Trapped Fermi gases is rapidly emerging as a new important area of research in the field of ultracold atomic gases. 
Impressive experimental results have been presented with the demonstration of trapping and cooling of fermionic 
alkalis using a magnetic trap for 40 K Q and 6 Li Q), and an optical trap for 6 Li These systems provide a large 
degree of experimental control: One can manipulate the effective particle-particle interaction by tuning an external 
magnetic field the geometry of the system can be changed, the number of internal degrees of freedom of the 
, constituent particles is given by the number of trapped hyperfine states, and the fundamental quantum statistics can 
be altered by trapping alkali atoms of Bose, Fermi (or both) statistics. Ultracold atomic Fermi gases therefore provide 
^sO , a very rich system to study with useful concepts coming from atomic, condensed matter, nuclear, and statistical 
physics and several theoretical papers have already been published focusing on various aspects of this subject 

The high precision spectroscopy of collective modes possible in these systems combined with the ability to ma- 
nipulate the interaction strength makes the study of collective modes for trapped Fermi gases a promising subject. 
Indeed, such studies have proven to be very useful in the field of Bose-Einstein condensation ||. The collective modes 
which can be observed as shape oscillations of the atomic cloud correspond to sound waves in a homogeneous gas, 
but for a confined system their spectrum is discrete. It is well known that there in general are two regimes in which 
to study these modes ||: When the lifetime t of the quasiparticles is much longer than the characteristic period 
I ■ of motion (i.e. lutt 3> 1 for atoms in a harmonic trap of frequency lot), there are few scattering events per sound 
oscillation, and the restoring force is due to the self-consistent mean field of the gas. Wave motion encountered in 
this limit is designated "zero sound". For the hydrodynamic regime ujtt «C 1, on the other hand, collisions ensure 
local thermodynamic equilibrium. Hydrodynamic sounds waves are sometimes called "first sound" . Due to the large 
variations in the effective interaction between the relevant alkali atoms in different hy per fine states, both regimes 
could be experimentally relevant and they have thus both been studied theoretically |l^,|ll| . 

We will in this paper examine the spectrum of collective modes in the collisionless regime for a gas in a spherical 
harmonic trap. By using a formalism well-known from nuclear physics, we systematically study the effect of the 
interactions on various aspects of the collective mode spectrum for modes of monopole, dipole, and quadrupole 
symmetry when the gas is in the normal phase. We show that for weak coupling, the spectrum is characterized 
by a large degree of collectivity (giant resonances). For such weak interactions where the shift of the collective 
mode spectrum away from the ideal gas result is very small, a sum rule calculation of the collective mode spectrum 
assuming a non-interacting ground state is shown to work well. With increasing coupling however, the effect of the 
interactions on the ground state become significant thus complicating a sum rule calculation. For stronger coupling, 
we demonstrate that effects such as Landau damping and mode splitting become important for some modes. The 
magnitude of these effects, which depends on the overlap between the single particle and the collective response, is 
shown to be different for the various modes considered. Finally, we briefly discuss the effects of a finite temperature 
on the spectrum, and higher modes relevant for the thermodynamics of the system are considered. 



II. FORMALISM 



In this section, we outline the basic theoretical framework used to calculate the collective mode spectrum of a 
trapped gas of interacting fermions. We consider a gas of fermionic atoms of mass m, confined by a potential Uq{t), 
with an equal number of atoms N a in each of two hyperfine states, \a =|, j). Two fermions in the same internal state 
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cr must have odd relative orbital angular momentum (minimally p-wave) , and at low temperatures T the centrifugal 
barrier suppresses their mutual interaction |1S^| . Thus, we assume the interaction to be effective only between atoms 
in different hyperfine states, and to be dominated by the low-energy s-wave scattering characterized by the scattering 
length a. Within mean field (Hartree Fock) theory, the gas is then described by the Hamiltonian: 

H = J2j rf 3 r^(r)[^V 2 + ^r 2 -H^(r)+3^ / d 3 r (^(r)^(r))4(r)^(r) + F(t), (1) 

(7 a 

where (...) denotes the thermal average, and the field operators obey the usual fermionic anticommutation rules 

describing the annihilation of a fermion at position r in the hyperfine state \a). The trapping potential is for simplicity 
assumed to be an isotropic harmonic oscillator Uo(r) — ^mw^r 2 , /ip is the chemical potential, and g — 4naH 2 /m. 
F{t) is an external perturbation with which we probe the system. The subject of this paper is the effects of the 
interaction on the collective mode spectrum of the gas in the normal phase and we are therefore neglecting pairing 
terms describing a possible superfluid transition. 

When the atom-atom interaction is rather weak such that the life-time r of the quasiparticles is much larger than 
1/Vr, the appropriate way to calculate the spectrum is the self-consistent random phase approximation (RPA) p3| . 
The problem of collective modes of a spatially confined two-component Fermi system has been studied extensively in 
the field of nuclear physics Q|. Thus, the basic formalism used in this paper is well-known from this field. 

As we are interested in modes corresponding to density fluctuations excited by perturbations of the kind F(t) oc 
exp(iwi) J2 a I d? r F<r( r )p<r( r )i we consider the retarded density-density correlation function {{p a (r)p a >(r')))(uj). Here, 
p a {r) = t/£ (r)-0o-(r) is the density operator for atoms in the hyperfine state \a) and ((AB))(uj) is the Fourier transform 
of the retarded function -i<d(t - t'){[A, B\) @. Within RPA, we have 

n(u) = [l-gB a (u)]- 1 %(u) (2) 

with 



fl(w) = 




fW^-<! «Mr)PT(r')»oM \ m 

n ° H ^ «p ; (r)p ; (r')»oM J (4) 
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and 

" ~ h \ S(r - r') 

Here the matrix products denote integrals: AB = J d 3 r"A(r, r")B(r", r'). The independent particle correlation 
functions ((. . . }}o are calculated within the time-independent Hartree-Fock (HF) approximation. In order to describe 
the effect of the interactions on the collective mode spectrum correctly, it is crucial that one uses the self-consistent 
static HF states to calculate the independent particle correlation functions ((. ..))o- If one simply uses the non- 
interacting trap states, the effect of the interactions on the single particle excitations are ignored and the predicted 
collective mode spectrum from Eq. (0) will be incorrect. 

For a spherically symmetric trap, the correlation functions split into terms describing the various multipole modes. 
Indeed, we have 

((pAv)p^{v')))(uj) - ^ W (ro ) rV» i y i „(M)F i M(fi^T (6) 

LM 

with r denoting the radial distance to the center of the trap and the Yz,m(0> 4 1 ) being the usual spherical harmonics. 
For the independent particle correlation functions, we obtain for the radial part p6[ : 

pp (ra,r a ,uj) l = > — (// 00 \LU)\ (J v l - J v n> ) ■ (") 

£rf v Att(2L + iy r ' 2 Z V 1 - Cfi' - »w - *a 

The quasiparticle wave functions u v i(r) are obtained from a self-consistent solution to mean field Hamiltonian in 
Eq.(l) with F(t) = writing the HF wave functions on the form u v i m (r) = u^\Xlm^i<^)It [9. The quasiparticles 
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u v im(^) have energies and f v i = [exp(/3£„/) + is the usual Fermi function with f3 = 1/ksT. We write (ll'00\L0) 
for the Clebsch-Gordan coefficients. 

The structure of the self-consistent RPA calculation is then the usual one: First we obtain a self-consistent solution 
to the mean field Hamiltonian in Eq. (1) with F = 0; then fto(cj) is formed from Eq.([|),(6)-(7) for varying external 
frequency oj; and finally the RPA response function is calculated from Eq. (||). The poles of Tl(u>) then give the 
collective modes of the gas excited by the density operator. 

It follows from Eq. (6) that the spherical symmetry permits an independent calculation of the radial part of the 
response functions for each spherical multipole, and the numerical dimension of the matrices is therefore determined 
by the number of radial mesh points. Hence, the method outlined is computationally very efficient as compared to a 
usual particle-hole representation of the RPA which in general yields matrices of very large dimensions jfffl . 



III. OBSERVABLES AND SUM RULES 



A quantity of great physical interest is the strength function directly related to the net transitions per unit time 
with energy hu induced by the operator F: 



e -0E n _ e -m 



S(F,u) = % \{n\F\m)\ 2 8(hw + E n - E m ). (8) 

nm 

Here Z is the grand partition function and \n) is an eigenstate of the Hamiltonian with energy E n . For operators of 
the form F(t) oc I d 3 rF a (r)p a (r), it can be obtained as 



S(F,cj) = -^J2 J d 3 rd 3 r'F a (r)F a/ (r')Im[n a ^(r,r\ 



(9) 



The strength function is the main object of the present paper. We will calculate it for perturbing operators of the 
form F a (r) = F a (r)YLM{d,Tr) exciting the monopole (L = 0), dipole (L — 1), and quadrupole (L = 2) modes. We 
take F| (r) = F± (r) for the monopole and quadrupole modes thereby exciting total density fluctuations of the given 
symmetries. For the dipole symmetry, we take F^(r) = ±i^(r) exciting both the dipole, and the spin-dipole mode. 
For the dipole mode, we have u> = lot independent of the interactions since it simply corresponds to a center-of-mass 
oscillation of the cloud. We calculate it to check the consistency of our numerics. 

If one defines the moments = J duiS(F, uj)ui k , the frequency u) of the mode excited by F can be calculated as 



« = (10) 
V rrik 

if the operator only excites one mode (giant resonance). This method is very handy since the two moments can be 
calculated as mi = (0| [F, [H, F}] |0) and — (0\\\F, H], \H, \H, F]]] |0); i.e. one only needs information about the 



ground state of the system to calculate the excitation energy |18 |. 

An important rule coming from particle conservation is the /-sum rule p6| : 



/,/, = J dwS(F,u)cj = y"d 3 r|VF CT (r)| 2 p CT (r) I LI ) 



with Po-(r) = (p a (r)) being the equilibrium density. This identity provides another important check on the reliability 
of our RPA calculations as the equilibrium density is easily obtained from the static HF solution. For our numerical 
RPA results, this sum rule turns out to be obeyed to within ~ 3% indicating the accuracy of our calculations. 



IV. RESULTS 



In this section, we present some typical results concerning the effect of the interactions on the collective mode 
spectrum obtained by solving the RPA equations for varying coupling strengths. 

As outlined above, if the perturbing operator only excites one state, then a sum-rule calculation based on Eq. jlfj| ) 
should give the correct energy of the excited state if an accurate ground state is used. Based on this assumption and 
on the fact that for weak interactions the ground state can be well approximated by an ideal Fermi gas, the lowest 
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collective mode of monopole, (spin)-dipole, and quadrupole symmetry was calculated in Ref. jTo|. Several effects can 
make such a calculation invalid: The interactions can change the ground state significantly away from the ideal gas 
limit although the operator F still excites a single collective mode (giant resonance). For stronger couplings, the 
mode can be fragmented into several modes and/or an incoherent excitation of several single particle states. This 
latter effect is known as Landau damping. We will now consider the importance of these effects for varying coupling 
strengths. 

For very weak interactions, we expect the sum-rule approach to give an accurate description of the lowest collective 
modes. We have checked this by comparing our RPA results for small coupling strengths with the prediction based 
on a sum rule calculation. For simplicity, we here only show results for the monopole mode [F a (r) = F(r)Yoo(6, (f))]. 
Completely equivalent results are obtained for the spin-dipole and the quadrupole modes. Using a local density 
approximation for the T = ground state of an ideal gas, a sum rule calculation predicts the frequency of the lowest 
monopole mode to be wm — 2u>t(1 + 0.3 x ^kpa) 1 ^ 2 with kp = y^2mfj,p/h 2 being the Fermi vector in the center 
of the trap |10|. Fig ure ([l|) d epicts a typical result obtained from a RPA calculation. We plot S(F,u>) (in units of 
l\/hujT with lh = y/h/mujT being the trap length) given by Eq. (g) for F a {r) — r 2 . The parameters are chosen 
such that there are ~ 10 4 atoms trapped and T — 0. Also, we have added a small imaginary part T — O.Q5hu>T 
to the frequency to model a smooth response typically observed experimentally The strength of the interaction is 
characterized by the parameter kpa. The response is calculated as a function of the frequency ui for the coupling 
strengths: kp\a\ ~ 0.06,0.13,0.19,0.25 and 0.32. We have chosen an attractive interaction (a < 0) such that the 
frequency shift away from the ideal gas result lum = 2lut is increasingly negative with increasing \a\. The prediction 
for lom based on the sum-rule approach is indicated by the symbol x for the various coupling strengths. 

As can be seen from Fig. (|IJ), there is good agreement between the prediction for the frequency of the mode based 
on the sum-rule approach and the RPA result for kp\a\ < 0.2. For such weak coupling strengths, the r 2 perturbation 
only excites the lowest monopole mode and the ground state is essentially given by the ideal gas result. For stronger 
interactions, the sum-rule result underestimates the frequency shift of the monopole mode due to the interactions. 
As can be seen from the inset in fig. (Q), the r 2 still essentially only excites one mode (giant resonance) for these 
coupling strengths. This means that a sum-rule calculation in principle should yield the correct result if one used to 
the exact ground state. We conclude that the discrepancy in fig. (|l|) is due to the fact that the simple non-interacting 
approximation to the ground state is inadequate in a sum-rule calculation for these coupling strengths. 

We have performed calculations with varying number of particles trapped and different coupling strengths for both 
the monopole, spin-dipole, and the quadrupole modes and compared with the predictions based on the sum-rule 
approach as given in Ref. |10[ . The general conclusion is, that the simple sum- rule results work well for kp\a\ < 0.2 
where the frequency shift away from the ideal gas result is very small for all these modes. For stronger interactions, one 
still has a giant resonance for the lowest monopole, spin-dipole, and quadrupole modes. However, to get an accurate 
prediction for the eigenfrequencies for these coupling strengths based on a sum-rule calculation, one would need to 
take into account the effect of the interactions on the ground state. Such a calculation could involve determining the 
ground state using a static Hartree Fock theory and then calculating the moments m^. 

We now consider the collective mode spectrum for more strongly interacting systems where the eigenmodes differ 
significantly from the ideal gas result. In this regime, the results based on the simple sum-rule calculation are invalid. 
This limit is especially relevant for experiments on 6 Li, where large an negative scatteringlengths between certain 
hyperfine states have been predicted theoretically and recently measured experimentally |3j. I n n g- (H) - ©! we P^°* 
the response S(F,lu) for the monopole, spin-dipole, and quadrupole symmetries excited by F a (r) — r 2 , F„ = —F s = r, 
and F a (r) = r 2 respectively. We also plot the dipole response excited by F a = v to further establish the accuracy 
of the calculations. There are ~ 1.6 x 10 4 particles trapped, kpa ~ —0.6, T — O.Obhui, and T — 0. The symbols x 



again indicate the sum-rule predictions for the collective modes |10 . We also plot the single particle response Sq(F,oj) 
based on using Ho(lu) instead of H(ui) in Eq.( ^). The collective response is indicated by solid lines whereas the single 
particle response is indicated by dotted lines. 

First, we note from fig. (||) that the dipole response is very accurately peaked at u> = ojp as it should be. This 
provides another confirmation of the accuracy of our calculations. Also, we see that the single particle response 
(dashed lines) in general is moved to higher frequencies as compared to the ideal gas result for all modes. This is 
because the mean field is attractive for a < and concentrated in the center of the trap. Since the spatial extend of 
the trap states in general increases with their energy, the decrease in the corresponding eigen-energies due to the mean 
field is larger for the lower trap states compared to for the higher states. The gap between the harmonic oscillator 
bands thus increases compared to the non-interacting limit resulting in a shift of the single particle response to higher 
frequencies [p||l9|]. Correspondingly, for a repulsive interaction the shift in the single particle response would have 
been to lower frequencies. We see that for all modes (except, of course, for the dipole mode), there is a significant 
disagreement between the sum-rule results and the RPA prediction as expected. 

Figure (0) indicates that the monopole mode is still characterized by a single giant resonance excited by the operator 
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r 2 . This is because contrary to the spin-dipole and quadrupole modes, the collective response for the monopole mode is 
moved to lower frequencies compared to the ideal gas result (u>m = 2u>t) due to the interactions. As can be seen from 
fig. (||), the overlap between the collective response and the single particle response (which is moved to frequencies 
lu > 2lut) is therefore insignificant prohibiting a decay of the collective mode into single particle excitations. The 
monopole mode is therefore well-defined. 

The strength of spin-dipole mode depicted in fig. (|^) is now distributed over a large frequency region. Both the 
collective and the single particle response is shifted to lu > lut resulting in a significant overlap between the collective 
mode and the single particle excitations. The collective mode can therefore decay into an incoherent motion of 
individual quasi particles (Landau damping). We also note from fig. (||) that this overlap with the single particle 
spectrum gives rise to a second small resonance located at lu ~ 1.135ojt- Obviously, such behaviour cannot be 
described within a sum-rule approach. The general shift to a higher frequency for the spin-dipole collective response 
as compared to the ideal case (w = lot) directly follows from the fact that for a < 0, the mean field provides an extra 
attractive force for the two hyperfine species oscillating in anti-phase and the spring constant for this motion is thus 
enhanced. 

For the quadrupole mode, there is a very low frequency response in the region < u> ^ 0.2lot as is highlighted in the 
inset in fig. (^|). This low frequency response is straightforward to understand: The quadrupole mode is characterized 
by single particle excitations with angular momentum change Al = ±2. The non-interacting spherical harmonic 
oscillator energies are given by the bands E n i = (n + 3/2)tuo with I = 0, 2, . . .n for n even and I = 1, 3, . . . n — 1 for n 
odd. Hence, the low frequency part of the quadrupole response is simply due to quasi particle excitations within such 
an energy band. If the gas was ideal, the frequency of such mira-band transitions would be and there would be no 
net transitions. Due to mean field effects there is a slight dispersion within each band |(| resulting in a finite response 
at small lo. For the monopole and spin-dipole modes, we have Al = and Al = ±1 and there are no intra-band 
transitions excited. These two modes therefore do not have any low frequency response until the interactions are so 
strong as to make each harmonic oscillator band overlap such that inter-band transitions have a low frequency. This 
low frequency response should have a significant effect on the quadrupole compressibility of the gas since from simple 
1. order perturbation theory one easily sees that the effect of a static quadrupole perturbation on the shape of the 
cloud is determined by the low energy excitations. 

The gas is in principle supcrfluid for T = for an attractive interaction, and fcsT c ~ 2.81ilot for the parameters 
given above Q. The presence of superfluidity could change significantly the low frequency response of the gas and 
this effect will be examined in a future publication. Since the gas strictly is superfluid for T — for an attractive 
interaction, we need for consistency to examine whether the effects described above relevant for a gas in the normal 
phase are still present for T > T c . In Fig. (||), we therefore also plot the collective response for ksT — 3thoT > ksTc 
indicated as a dotted line. We see that the effect of a finite temperature on the spectrum is as expected to smear 
the response out somewhat due to the thermal occupation of excited states. Furthermore, as a finite T tends to 
make the gas more dilute and thus diminish the effect of the interactions, the spectrum is shifted slightly towards the 
ideal gas limit. However, the effect is quite small for the low T considered here and the main features of the T = 
spectrum remain. Importantly, we note that even though UbT — 3Hlut, the low frequency response (lo < 0.2lut) is 
not completely suppressed. Obviously, for increasing T the response the gas eventually becomes ideal. 

In order to excite higher modes of the gas, one in general needs perturbations containing higher powers of r as can 
be seen if one writes these operators in terms of the raising and lowering operators for the harmonic oscillator. This 
is illustrated in fig. (|j), where we plot the response to a perturbation given by F a (r) — r 4 , F a = —F s = r 3 , and 
F a (r) = r 4 for the monopole, spin-dipole and quadrupole respectively. We have chosen ~ 1.6 x 10 4 particles trapped, 
kpa ~ —0.6, T = O.lhw, and T — 0. We see that for such perturbations, higher modes are also excited as expected. 
However, the strength of these higher excitations is rather low compared to the lowest modes. This is due to the 
fact that these modes correspond to additional nodes in the density fluctuation profiles and their overlap with the 
perturbing operator F a (r) is therefore smaller than for the lowest modes. In order to increase the excitation rate of 
these higher modes, one could introduce a perturbation with more such radial nodes. In general, we believe that the 
excitation of these higher modes is somewhat more complicated than the excitation of the lowest modes, as one needs 
to experimentally generate spatial perturbations of a non-harmonic kind. 

V. CONCLUSION 

We have systematically examined the effects of the particle-particle interactions on the collective mode spectrum 
of a normal trapped Fermi gas within a self-consistent RPA scheme. For weak interactions where the shift in the 
eigenfrequencies away from the ideal case is very small, it was shown that a simple sum-rule approach to calculate 
the frequencies of the lowest modes works well. With increasing interaction strength however, such a calculation was 
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demonstrated to be complicated by the fact that it is necessary to include the effects of the particle-particle scattering 
on the ground state to achieve accurate results. As the coupling strengths increases further, we demonstrated the 
emergence of effects such as mode splitting and Landau damping coming from the interplay between the collective and 
the single particle excitations. The importance of such effects was shown to depend on the symmetry of the mode. 
For the quadrupole mode, a very low frequency response was predicted with possible implications on the quadrupole 
compressibility of the gas. The temperature dependence of the modes was considered and we finally briefly discussed 
the excitation of some of the higher collective modes. Due to the possibility of experimentally manipulating the 
effective interaction between alkali atoms and the high precision spectroscopy typically obtainable in such systems, we 
believe that some of the effects described in this paper should be observable using present day experimental technology. 
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FIG. 1. The T = collective response S(F,u>) in units of l^/hur for the monopole mode excited by F(r) = r 2 for various 
coupling strengths. The sum-rule results are connected with the corresponding RPA results by an arrow. The inset shows 
S(F,ui) for fci?H — 0.32 over a larger frequency range. 
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FIG. 2. The T = collective (solid line) and single particle (dashed line) response for the monopole (in units of I^/Hlot) and 
dipole mode (in units of I^/Hojt) excited by F a {r) — r 2 and F a (r) — r respectively. The coupling strength is fcjr|o| ~ 0.6. 
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Bruun figure 3 
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FIG. 3. The T — collective (solid line) and single particle (dashed line) response for the spin-dipole (in units of l\/hujT) 
and the quadrupole mode (in units of l^/hiOr) excited by F a (r) = — Fair) = r and F a (r) = r 2 respectively. The coupling 
strength is fe_p|a| — 0.6. For the quadrupole mode, the dotted line depicts the collective response for ksT = 3hoJT- The inset 
shows in more detail the low frequency collective quadrupole response. 




FIG. 4. The T = collective (solid line) and single particle (dashed line) response for the monopole (in units of I^/Hljt), 
spin-dipole (in units of 1%/Tujt), and quadrupole modes (in units of I^/Hujt) excited by F a (r) = r 4 , F a (r) = —Fair) — r 3 and 
F a (r) = r 4 respectively. 
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